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1. Preamble 

Hydrogen in the universe recombined about half a million years after the Big Bang, and 
cooled down to a temperature of a few kelvins until the first non-linearities developed, and 
evolved into stars, galaxies, and black holes that lit up the universe again. In currently 
popular cold dark matter flat cosmologies (ACDM), some time beyond a redshift of 10 
the gas within halos with virial temperatures Tvir ^ 10"* K - or, equivalently, with masses 
M > 10^ [(l + 2;)/10]~'^/^ Mq - cooled rapidly due to the excitation of hydrogen Lya and 
fragmented. Massive stars formed with some initial mass function (IMF), synthesized 
heavy elements, and exploded as Type II supernovae after a few xlO^yr, enriching 
the surrounding medium: these subgalactic stellar systems, aided perhaps by an early 
population of accreting black holes in their nuclei, generated the ultraviolet radiation 
and mechanical energy that contributed to the reheating and reionization of the cosmos. 
It is widely believed that coUisional excitation of molecular hydrogen may have allowed 
gas in even smaller systems - virial temperatures of a thousand K, corresponding to 
masses around 5 x 10^ [(1 + z)/10]~^/^ M0 - to cool and form stars at even earlier times 
(Couchman & Rees 1986; Haiman et al. 1996; Tegmark et al. 1997; Abel et al. 2000, 
2002; Fuller & Couchman 2000; Bromm et al. 2002; Reed et al. 2005; Kuhlen & Madau 
2005). Throughout the epoch of structure formation, the all-pervading intergalactic 
medium (IGM), which contains most of the ordinary baryonic material left over from 
the Big Bang, becomes clumpy under the influence of gravity, and acts as a source for 
the gas that gets accreted, cools, and forms stars within subgalactic fragments, and as a 
sink for the metal enriched material, energy, and radiation which they eject. The well- 
established existence of heavy elements like carbon and silicon in the Lya forest clouds 
at 2; = 2 — 6 (Songaila 2001; Pettini et al. 2003; Ryan- Weber et al. 2006) may be indirect 
evidence for such an early episode of pregalactic star formation. The recently released 
Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data require the universe to be 
fully reionized by redshift z^. = H ± 2.5 (Spergel et al. 2006), another indication that 
significant star-formation activity started at very early cosmic times. 

The last decade has witnessed great advances in our understanding of the high redshift 
universe, thanks to breakthroughs achieved with satellites, 8-lOm class telescopes, and 
cosmic microwave background (CMB) experiments. Large surveys such as the Sloan 
Digital Sky Survey (SDSS), together with the use of novel instruments and observational 
techniques have led to the discovery of galaxies and quasars at redshifts in excess of 6. At 
the time of writing, 9 quasars have already been found with z > 6 (Fan 2006), and one 
actively star- forming has been spectroscopically confirmed at z = 6.96 (lye et al. 2006). 
These sources probe an epoch when the universe was < 7% of its current age. Keck 
and VXT observations of redshifted H i Lya ('forest') absorption have been shown to be 
sensitive probe of the distribution of gaseous matter in the universe (see Ranch 1998 for 
a review). Gamma-ray bursts have recently displayed their potential to replace quasars 
as the preferred probe of early star formation and chemical enrichment: GRB050904, the 
most distant event known to date, is at z = 6.39 (Haislip et al. 2006). The underlying 
goal of all these efforts is to understand the growth of cosmic structures, the properties 
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of galaxies and their evolution, and ultimately to map the transition from the cosmic 
"dark age" to a ionized universe populated with luminous sources. 

Progress has been equally significant on the theoretical side. The key idea of currently 
popular cosmological scenarios, that primordial density fluctuations grow by gravitational 
instability driven by coUisionless CDM, has been elaborated upon and explored in detail 
through large-scale numerical simulations on supercomputers, leading to a hierarchical 
('bottom-up') scenario of structure formation. In this model, the flrst objects to form are 
on subgalactic scales, and merge to make progressively bigger structures ('hierarchical 
clustering'). Ordinary matter in the universe follows the dynamics dictated by the dark 
matter until radiative, hydrodynamic, and star formation processes take over. According 
to these calculations, a truely inter- and proto-galactic medium (the main repository 
of baryons at high redshift) collapses under the influence of dark matter gravity into 
flattened and filamentary structures, which are seen in absorption against background 
QSOs. Gas condensation in the first baryonic objects is possible through the formation 
of H2 molecules, which cool via roto-vibrational transitions down to temperatures of a 
few hundred kelvins. In the absence of a UV photodissociating fiux and of ionizing X- 
ray radiation, three-dimensional simulations of early structure formation show that the 
fraction of cold, dense gas available for accretion onto seed black holes or star formation 
exceeds 20% for halos more massive than 10^ Mq already at redshifts 20 (Machacek et al. 
2003; Yoshida et al. 2003). 

In spite of some significant achievements in our understanding of the formation of 
cosmic structures, there are still many challenges facing hierarchical clustering theories. 
While quite successful in matching the observed large-scale density distribution (like, e.g., 
the properties of galaxy clusters, galaxy clustering, and the statistics of the Lya forest), 
CDM simulations appear to produce halos that are too centrally concentrated compared 
to the mass distribution inferred from the rotation curves of (dark matter-dominated) 
dwarf galaxies, and to predict too many dark matter subhalos compared to the number of 
dwarf satellites observed within the Local Group (Moore et al. 1999; Klypin et al. 1999). 
Another perceived difficulty (arguably connected with the "missing satellites problem" , 
e.g. Bullock et al. 2000) is our inability to predict when and how the universe was 
reheated and rcionizcd, i.e. to understand the initial conditions of the galaxy formation 
process and the basic building blocks of today's massive baryonic structures. We know 
that at least some galaxies and quasars had already formed when the universe was less 
than 10^ yr old. But when did the first luminous clumps form, was star formation 
efficient in baryonic objects below the atomic cooling mass, and what was the impact 
of these early systems on the surrounding intergalactic gas? The crucial processes of 
star formation and "feedback" (e.g. the; efl:ect of the energy input from the earliest 
generations of sources on later ones) in the nuclei of galaxies are still poorly understood. 
Accreting black holes can release large amounts of energy to their surroundings, and 
may play a role in regulating the thermodynamics of the interstellar, intracluster, and 
intergalactic medium. The detailed astrophysics of these processes is, however, unclear. 
Although we may have a sketchy history of the production of the chemical elements in the 
universe, we know little about how and where exactly they were produced and how they 
are distributed in the IGM and in the intracluster medium. Finally, where are the first 
stars and their remnants now, and why are the hundreds of massive satellites predicted 
to survive today in the Milky Way halo dark? 

In these lectures I will describe some of the basic tools that have been developed to 
study the dawn of galaxies, the lessons learned, and summarize our current understanding 
of the birth of the earliest astrophysical objects. 
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2. The Dark Ages 

2.1. Cosmological preliminaries 

Recent CMB experiments, in conjunction with new measurements of the large-scale struc- 
ture in the present-day universe, the SNIa Hubble diagram, and other observations, 
have lead to a substantial reduction in the uncertanties in the parameters describing 
the ACDM concordance cosmology. Wc appear to be living in a flat {il„i + iljy = 1) 
universe dominated by a cosmological constant and seeded with an approximately scale- 
invariant primordial spectrum of Gaussian density fluctuations. A ACDM cosmology 
with 0,„ = 0.24, ilA = 0.76, fli, = 0.042, h = i/o/lOO kms-^Mpc-^ = 0.73, n = 0.95, 
and o"8 = 0.75 is consistent with the best-fit parameters from the WMAP 3-year data 
release (Spergel ct al. 2006). Here fl„i = Pm/Pc^ the present-day matter density 
(including cold dark matter as well as a contribution from baryons) relative to the 
critical density p° = SHq/SttG, fl^ is the vacuum energy contribution, Hq is the Hubble 
constant today, n is the spectral index of the matter power spectrum at inflation, and trg 
normalizes the power spectrum: it is the root-mean-square amplitude of mass fluctua- 
tions in a sphere of radius 8h~^ Mpc. The lower values of as and n compared to WMAP 
1-year results (Spergel et al. 2003) have the effect of delaying structure formation and 
reducing small-scale power. In this cosmology, and at the epochs of interest here, the 
expansion rate H evolves according to the Friedmann equation 

H = ~ = Ho[n^{i + zf + nA]'/^ (2.1) 

where z is the redshift. Light emitted by a source at time t is observed today (at time 
to) with a redshift z = l/a^t) — 1, where a is the cosmic scale factor [a{to) = 1]. The age 
of the (flat) universe today is 



*° " ' (l + z')H{z') ~ S^oTOI 



1.025 Fq-^ = 13.7 Gyr (2.2) 



(Ryden 2003). At high redshift, the universe approaches the Einstein-de Sitter behaviour, 
a oc f^/^, and its age is given by 

t{z)^ ?=(l + ^)-3/2. (2.3) 

3i/ov "to 

The average baryon density today is ri° = 2.5 x 10~^ cm~^, and the hydrogen density 
is riH = (1 — Yp) Ub = 0.75 where Yp is the primordial mass fraction of helium. The 
photon density is n° = {2A/tt'^) {kBTa/hc)'^ = 400 cm'^, where Tq = 2.728 ± 0.004 
K (Fixsen et al. 1996), and the cosmic baryon-to- photon ratio is then t] = nb/uj = 
6.5 X 10~^°. The electron scattering optical depth to redshift is given by 

re{zr)=j^ dzn^arc- = J^ (TT^)^ ^'"'^ 

where Ug is the electron density at redshift z and ax the Thomson cross-section. Assum- 
ing a constant electron fraction Xe = ne/nn with redshift, and neglecting the vacuum 
energy contribution in H{z), this can be rewritten as 

Te{Zr) = ""f^^ r dzVT+^ = 0.0022 [(1 + Zrf/^ - 1]. (2.5) 
J^oV^'m Jo 



t Hereafter densities measured at the present epoch will be indicated either by the superscript 
or the subscript '0'. 
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Ignoring helium, the 3-year WMAP polarization data requiring Te = 0.09 ± 0.03 are 
consistent with a universe in which Xe changes from essentially close to zero to unity at 
Zr — Hi: 2.5, and Xe ^ 1 thereafter. 

Given a population of objects with proper number density n{z) and geometric cross- 
section sigma the incremental probablity dP that a line of sight will intersect one 
of the objects in the redshift interval dz at redshift z is 

dP = n[zmz)cf^dz = nizmzy^^^y (2.6) 

2.2. Physics of recombination 

Recombination marks the end of the plasma era, and the beginning of the era of neutral 
matter. Atomic hydrogen has an ionization potential of 7 = 13.6 eV: it takes 10.2 eV to 
raise an electron from the ground state to the first excited state, from which a further 
3.4 eV will free it, 

7 + H< — >p + e-. (2.7) 

According to Boltzmann statistics, the number density distribution of non-relativistic 
particles of mass rrii in thermal equilbrium is given by 

, rp\ 3/2 / 2 



where Qi and jjii are the statistical weight and chemical potential of the species. When 
photoionization equibrium also holds, then /Xg + = Mh- Recalling that = gp = 0.5 
and gu = 2, one then obtains the Saha equation: 

neUp ^ ^ {2nmekBTf/^ (-^] , (2.9) 

nnnm I - Xe nn{27rh)^ \ ksT J ' 

which can be rewritten as 

^.2 



In ^ = 52.4 - 1.5 ln(l + z) - 58, 000 (1 + z)-\ (2.10) 

\i--XeJ 

The ionization fraction goes from 0.91 to 0.005 as the redshift decreases from 1500 to 
1100, and the temperature drops from 4100 K to 3000 K. The elapsed time is less than 
200,000 yr. 

Altliongli the Saha equation describes reasonably well the initial phases of the depar- 
ture from complete ionization, the assumption of equilibrium rapidly ceases to be valid 
in an expanding universe, and recombination freezes out. The residual electron fraction 
can be estimated as follows. The rate at which electrons recombine with protons is 

dne _ np 

— = -aBUeUp = -- — , (2.11) 
dZ ^rec 

where irec is the characteristic time for recombination and as is the radiative recom- 
bination coefficient. This is the product of the electron capture cross-section cr„ and 
the electron velocity Ve, averaged over a thermal distribution and summed over all ex- 
cited states n > 2 of the hydrogen atom, as = '^{(JnVe)- The radiative recombination 
coefficient is well approximated by the fitting formula 

as = 6.8 X 10"^^ Tg-^-^ cm^ s"^ = 1.85 x 10"^° (1 + z)-°-^ cm^ s'S (2.12) 

where T3 = T/SOOO K and in the second equality I used T = ro(l + z). When the 
recombination rate falls below the expansion rate, i.e. when tree > ^/H, the formation of 
neutral atoms ceases and the remaining electrons and protons have negligible probability 
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Figure 1. Helium and hydrogen recombination for a WMAP 3-year cosmology. The step at 
earlier times in the left panel is due to the recombination of He iii into He ii. We used the 
code RECFAST (Seager et al. 1999) to compute the electron fraction Xe- Note that the residual 
electron fraction (determined by the condition trccH = 1) scales with cosmological parameters 
as Xe oc y/Tl^/flth. 



for combining with each other: 



XeU^aB 



, 0.0335, , 7 2.5x10-4 ^^^^^ 
(l + z)-3/2^ (l + 2)-o-7 = , (2.13) 



where I assumed z — 1100 in the third equality. The recombination of hydrogen in 
an expanding universe is actually delayed by a number of subtle effects that are not 
taken into account in the above formulation. An captured to the ground state of 
atomic hydrogen produces a photon that immediately ionizes another atom, leaving no 
net change. When an is captured instead to an excited state, the allowed decay to 
the ground state produces a resonant Lyman series photon. These photons have large 
capture cross-sections and put atoms in a high energy state that is easily photoionized 
again, thereby annulling the effect. That leaves two main routes to the production of 
atomic hydrogen: 1) two-photon decay from the 2s level to the ground state; and (2) 
loss of Lya resonance photons by the cosmological redshift. The resulting recombination 
history was derived by Peebles (1968) and Zel'dovich et al. (1969). 

In the redshift range 800 < z < 1200, the fractional ionization varies rapidly and is 
given approximately by the fitting formula 

/ Z \ 11-25 

= ( looo) ■ ^'-''^ 

This is a fit to a numerical output from the code RECFAST - an improved calculation 
of the recombination of H i. He i, and He ii in the early universe involving a line-by-line 
treatment of each atomic level (Seager et al. 1999). Using this expression, we can again 
compute the optical depth of the universe for Thomson scattering by free electrons: 

^ ' Jo {l + z)H{z) VlQQOy ^ ^ 

This optical depth is unity at Zdcc = 1100. From the optical depth we can compute the 
visibility function P{t), the probability that a photon was last scattered in the interval 
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Figure 2. Visibility function. Solid Ime: RECFAST with WMAP 3-year cosmological 
parameters. Dashed line: Jones & Wyse (1985). 



(z, z + dz). This is given by 

P(r) = ^ = 0.0038 (y^)" '' e-«-3(-/i°«°)^^ ^^ (2.16) 

and has a sharp maximum at z = Zdcc and a width of Az w 100. The finite thickness of the 
last scattering surface has important observational consequences for CMB temperature 
fluctuations on small scales. 

2.3. Coupling of gas and radiation 

The CMB plays a key role in the early evolution of structures. First, it sets the epoch of 
decoupling, when baryonic matter becomes free to move through the radiation field to 
form the first generation of gravitationally bound systems. Second, it fixes the matter 
temperature that in turns determines the Jeans scale of the minimum size of the first 
bound objects. 

Following Peebles (1993), consider an electron e~ moving at non-relativistic speed 
V <^ c through the CMB. In the e~ rest-frame, the CMB temperature measured at angle 
6 from the direction of motion is 

T{9) =T (^l + -cos9y (2.17) 
The radiation energy density per unit volume per unit solid angle dfl — d(j)d cos 9 around 

e is 

du=\aBT^i9)]— (2.18) 
An 

(here as is the radiation constant), and the net drag force (i.e. the component of the 
momentum transfer along the direction of motion, integrated over all directions of the 
radiation) felt by the electron is 

F = J arducose^ J ariasT^) [l + -cos 9 j cos9— 

4 arasT'^ 



3 c 



V. (2.19) 
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This force will be communicated to the protons through electrostatic coupling. The 
formation of the first gravitationally bound systems is limited by radiation drag, as the 
drag force per baryon is XgF, where Xe is the fractional ionization. The mean force 
divided by the mass nip of a hydrogen atom gives the deceleration time of the streaming 
motion: 

^_,^iyv^4aranT,Hl + z)\^ (2.20) 
V at 3 nipC 
The product of the expansion rate H and the velocity dissipation time tg is 

tsH = 7.6 X lO'^ hx-\l + z)-^/'^. (2.21) 

Prior to decoupling tsH <C 1. Since the characteristic time for the gravitational growth 
of mass density fluctuations is of the order of the expansion timescale, baryonic density 
fluctuations become free to grow only after decoupling. 

We can use the above results to find the rate of relaxation of the matter temperature Te 
to that of the radiation. The mean energy per electron in the plasma is i? = 3kBTp/2 = 
meiv"^)- The rate at which an electron is doing work against the radiation drag force if 
Fv, so the plasma transfers energy to the radiation at the mean rate, per electron: 

At thermal equilibrium Te = T, and this rate must be balanced by the rate at which 
photons scattering off electrons increase the matter energy: 

dE _4 araBT^ Sfcs ,rp rp^ ,^ 

Thus the rate of change of the matter temperature is 

1 — = ^ ^ ^ (T — T ). (2.24) 

dt iks dt (1 + Xe) 3meC ^ 

The factor Xe/{1 + Xe) accounts for the fact that the plasma energy loss rate per unit 
volume is —UedE/dt = —XeiindE/dt, while the total plasma energy density is {ne + np + 
nm)3kBTe/2 = nu{i + Xe)3kBTe/2. The expression above can be rewritten as 

dTe ^ T-Te 

dt tc ' 
where the "Compton cooling timescale" is 



(2.25) 



_ 3meC 1 + Xe _7Ax 10^^ s / 1 + a,e , 
' " AaraBT^ " (1 + zY \~2^ ' " ^ ' 



The characteristic condition for thermal coupling is then 



-1 



For fully ionized gas = 1 and the Compton cooling timescale is shorter than the 

expansion time at all rcdshifts z > Zc = 5- With increasing rcdshift above 5, it becomes 
increasingly difficult to keep optically thin ionized plasma hotter than the CMB. At 
redshift z > 100, well before the first energy sources (stars, accreting black holes) turn 
on, hydrogen will have the residual ionization Xg — 2.5 X 10 With this ionization, the 
characteristic relation for the relaxation of the matter temperature is 
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The coefficient of tlie fractional temperature difference readies unity at the "thermahza- 
tion redshift" zth ~ 130. That is, the residual ionization is enough to keep the matter in 
temperature equilibrium with the CMB well after decoupling. At redshift lower than zth 
the temperature of intergalactic gas falls adiabatically faster than that of the radiation, 
Tf, (X a^^. From the analysis above, the rate of change of the radiation energy density 
due to Compton scattering can be written as 



du 4 aTdBT^ Sfc^rif, 
dt Z c rrie 



(Te - T), (2.29) 



or 



— — My, dy = [rieCrTcdt) ^ — dr^ ^ . (2.30) 

Compton scattering causes a distorsion of the CMB spectrum, depopulating the Rayleigh- 
Jeans regime in favor of photons in the Wien tail. The "Compton-parameter" 

y = / j-^c?^ 2.31) 

is a dimensionless measure of this distorsion, and is proportional to the pressure of the 
electron gas rieksTe. The COBE satellite has shown the CMB to be thermal to high 
accuracy, setting a limit y < 1.5 x 10~^ (Fixsen et al. 1996). This can be shown to imply 

(xeTe) [(1 + zf^^ - 1] < 4 X 10^ K. (2.32) 

A universe that was reionized and reheated at z = 20 to {xe,Te) = (1, > 4 x 10^ K), for 
example, would violate the COBE y-limit. 
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reaction rate (cm"^ s ^ or s ^) notes 



H+ + e ^ H + 7 i?c2 = 8.76 x 10-"(1 + z)-°-^^ 

H + 7 ^ H++ e R2c = 2.41 x lO^^T^-^ exp (-^^) x i?c2 

H + e ^ H- + 7 1.4 X 10-i«T," '^28cxp (--^) 

H- + 7 ^ H + e 1.1 X 10-iT2 i3exp (-^) 

H- + H ^ H2 + e 1.5 X 10-9 Te < 300 

4.0 X IQ-^T-o i^ Te > 300 

H- + H+ ^ 2H 5.7 X lO-^T-o-s + 6.3 x IQ-^- 

9.2 X 10-"T0-5 + 4.4 X lO'^^Te 
H + H+ ^ H+ + 7 dcx[- 19.38 - 1.523 log Te+ 

1.118(logTe)2 - 0.1269(logre)3] < 10*-^ 

+ 7 ^ H + H+ 2.0 X lO^Ti-SS exp (-^fpo) v = 

1.63X 107exp(-2^) LTE 
H+ + H ^ H2 + H+ 6.4 X IQ-io 

H2 + H+ ^ H+ + H 3.0 X 10^10 exp (-^^) < 10^ 

1.5 X 10-1° exp (-™) Te > 10* 

He + H+ ^ HeH+ +7 7.6 x lO-^^T-O-S < 10^, 

3.45 X 10-16T-1-06 Te > 10^ 

HeH+ + H ^ He + H+ 9.1 x 10"^° 
HeH+ + 7 ^ He + H+ 6.8 x lO-^T^-^exp (-'^^) 



2.4. Hydrogen molecules in the early universe 

In the absence of heavy metals or dust grains, the processes of radiative cooling, cloud 
collapse, and star formation in the earliest astrophysical objects were undoubtedly quite 
different from those today. Saslaw & Zipoy (1967) and Peebles & Dicke (1968) were the 
first to realize the importance of gas phase reactions for the formation of the simplest 
molecule, H2, in the post-recombination epoch. The presence of even a trace abundance 
of H2 is of direct relevance for the cooling properties of primordial gas. In the absence of 
molecules this would be an extremely poor radiator: cooling by Lya photons is ineffective 
at temperatures < 8000 K, well above the matter and radiation temperature in the post- 
recombination era. It is the ability of dust-free gas to cool down to low temperatures 
that controls the formation of the first stars in subgalactic systems. 

The formation of H2 in the intergalactic medium is catalyzed by the residual free 
electrons and ions, through the reactions 

H + e- ^ R- + "f, (2.33) 
H- H ^ H2 + e-, (2.34) 

and 



H+ + H ^ H+ + 7, 
H+ + H ^ H2 + H+. 



(2.35) 
(2.36) 
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The second sequence produces molecular hydrogen at a higher redshift than the first, as 
the dissociation energy thresholds of the intermediaries species and H~ are 2.64 and 
0.75 eV, respectively. The direct radiative association of H + H ^ H2 + 7 is highly 
forbidden, owing to the negligible dipole moment of the homonuclear H2 molecule. The 
amount of H2 that forms depends on the evolution of the gas density. Galli & Palla (1998) 
have analyzed in detail the case where the hydrogen follows the general expansion of the 
universe. Table 1 shows the reaction rates for hydrogen and helium species that enter 
in what they call the minimal model, i.e. the reduced set of processes that reproduces 
with excellent accuracy the full chemistry of H2 molecules. Coltunn 1 gives the reaction, 
column 2 the rate coefficient (in cm^ s~^ for collisional processes, in s~^ for photo- 
processes), and column 3 the temperature range of validity of the rate coefficients and 
remarks on the rate (see Galli & Palla 1998). The gas and radiation temperatures are 
indicated by Te and T, respectively. 

The evolution of the abundance of H2 follows the well known behaviour (e.g. Lepp 
& Shull 1984) where an initial steep rise at z ~ 400 is determined by the channel, 
followed by a small contribution from H~ at 2; ~ 100. The freeze-out primordial fraction 
of H2 is [H2/H] ~ 10~^. This is too small to trigger runaway collapse, fragmentation, 
and star formation. The fate of a subgalactic system at high redshifts depends therefore 
on its ability to rapidly increase its H2 content following virialization. In the low-density 
limit (riH < 0.1 cm"'^), the H2 cooling function (in erg cm"' s^^) is well approximated 
by the expression 

logAnJnH ^ 0] = -103.0 -|- 97.59 log Te - 48.05(logTe)2 

+10.80(logTe)3 - 0.9032(logTe)^, (2.37) 
over the range 10 K < Te < 10* K (Galh & Palla 1998). 

3. The Emergence of Cosmic Structure 

3.1. Linear theory 

As shown above, it is only after hydrogen recombination that baryons can start falling 
into the already growing dark matter perturbations. Gas pressure can resist the force of 
gravity, and small-scale perturbations in the baryonic fluid will not grow in amplitude. 
However, at sufficiently large scales, gravity can overpower pressure gradients, thereby 
allowing the perturbation to grow. 

The linear evolution of sub-horizon density perturbations in the dark matter-baryon 
fluid is governed in the matter-dominated era by two second-order differential equations: 

5dm + 2H5i^ = \H^^l^ (/dm'Jdm + hh) (3-38) 

for the dark matter, and 

h + 2H6b = ^H^n'm (/dm^dm + fbSb) ' ^k^Sb (3.39) 

^ CI 

for the baryons, where 6dm{k) and db{k) are the Fourier components of the density 
fluctuations in the dark matter and baryons, f /dm and fb are the corresponding mass 
fractions, Cg is the gas sound speed, k the (comoving) wavenumber, and the derivates are 
taken with respect to cosmic time. Here = 8'KGp{t)/3H'^ = Qm{^+z)^/['^m{'^ + z)^ + 

f For each fluid component (i = b, dm) the real space fluctuation in the den- 
sity field, 5i(x) = Spi{x)/pi, can be written as a sum over Fourier modes, 
5i(x) = / d*k(27r)"^5i(k) expik-x. 
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(5dm = 47rGp5dm = (3.40) 



Qa] is the time-dependent matter density parameter, and p{t) is the total background 
matter density. Because there is 5 times more dark matter than baryons, it is the former 
that defines the pattern of gravitational wells in which structure formation occurs. In 
the case where fb — and the universe is static {H = 0), equation (3.38) becomes 

dyn 

where tdyn denotes the dynamical timescale. This equation admits solution Jdm = 

Ai cxp(f/fdyn) + ^2 cxp(— i/fdyn)- Aftcr a few dynamical times, only the exponentially 
growing term is significant: gravity tends to make small density fiuctuations in a static 
pressureless medium grow exponentially with time. The additional term oc HSdm present 
in an expanding universe can be thought as a "Hubble friction" term that acts to slow 
down the growth of density perturbations. Equation (3.38) admits the general solution 
for the growing mode: 

S,M = ^H^ohJ^^-^, (3.41) 

where the constant have been choosen so that an Einstein-de Sitter universe (fi^ = 1) 
Oa = 0) gives the familiar scaling (5dm (a) = a with coefficient unity. The right-hand side 
of equation (3.41) is called the linear growth factor D(a). Different values of ilm, lead 
to different linear growth factors: growing modes actually decrease in density, but not 
as fast as the average universe. Note how, in contrast to the exponential growth found 
in the static case, the growth of perturbations even in the case of an Einstein-de Sitter 
universe is just algebraic. A remarkable approximation formula to the growth factor in 
a flat universe follows from Lahav et al. (1991), 



<5dm(a) = D{a) 



^ 140 +140 '"^ 70 



(3.42) 



2(1 + z) 

This is good to a few percent in regions of plausible Qm, ^a- 

Equation (3.39) shows that, on large scales (i.e. small fc), pressure forces can be 
neglected and baryons simply track the dark matter fluctuations. On the contrary, on 
small scales (i.e. large fc), pressure dominates and baryon fluctuations will be suppressed 
relative to the dark matter. Gravity and pressure are equal at the characteristic Jeans 
scale: 

kj = —JAttGp. (3.43) 

Corresponding to this critical wavenumber kj there is a critical cosmological Jeans mass 
Mj, defined as the total (gas -|- dark matter) mass enclosed within the sphere of physical 
radius equal to na/kj, 

^^ = yH^J =yHi^G^J «8.8xl0^M,(-j , (3.44) 

where /i is the mean molecular weight. The evolution of Mj is shown in Figure 4. In the 
post-recombination universe, the baryon-electron gas is thermally coupled to the CMB, 
Tg oc a~^, and the Jeans mass is independent of redshift and comparable to the mass of 
globular clusters, Mj w IQ^ Mq. For z < Zth, the gas temperature drops as cx a~^, 
and the Jeans mass decreases with time, Mj oc a~"^/^. This trend is reversed by the 
reheating of the IGM. The energy released by the first collapsed objects drives the Jeans 
mass up to galaxy scales (Figure 4) : baryonic density perturbations stop growing as their 
mass drops below the new Jeans mass. In particular, photo-ionization by the ultraviolet 
radiation from the first stars and quasars would heat the IGM to temperatures of w 10^ K 
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(1+z) (1+z) 

Figure 4. Left: Same as Fig. 3 but for a universe reionized by ultraviolet radiation at 
Zr — 11. Right: Cosmological (gas + dark matter) Jeans mass. 



(corresponding to a Jeans mass Mj ~ 10^" M© at Zr ~ 11), suppressing gas infall into 
low mass systems and preventing new (dwarf) galaxies from forming. f 

3.2. Statistics of density fields 

The observed uniformity of the CMB guarantees that density fluctuations must have 
been quite small at decoupling, implying that the evolution of the density contrast can 
be studied at z < Zdoc using linear theory, and each mode S{k) evolves independently. 
The inflationary model predicts a scale-invariant primordial power spectrum of density 
fluctuations P{k) = (|i5(fc)P) oc fc", with n — 1 (the so-called Harrison- Zel'dovich spec- 
trum). It is the index n that governs the balance between large and small-scale power. 
In the case of a Gaussian random field with zero mean, the power spectrum contains 
the complete statistical information about the density inhomogeneity. It is often more 
convenient to use the dimensionless quantity Af. = [k^P{k)/2TT'^], which is the power 
per logarithmic interval in wavenumber k. In the matter-dominated epoch, this quantity 
retains its initial primordial shape (A| oc k"^^) only on very large scales. Small wave- 
length modes enter the horizon earlier on and their growth is suppressed more severely 
during the radiation-dominated epoch: on small scales the amplitude of is essentially 
suppressed by four powers of k (from /c"+'^ to fc""^). If n = 1, then small scales will have 
nearly the same power except for a weak, logarithmic dependence. Departures from the 
initially scale- free form are described by the transfer function T{k), defined such that 
T(0) = 1: 



P{k,z) = Afc" 



r^(fc), (3.45) 



D{0)_ 

where A is the normalization. An accurate fitting function for T(fc) in a ACDM universe 
is 

n = ^"^^+3^;^^^'^ [1 + 3.89g + (16.1?)^ + {^AQqf + (6.71?)^] , (3.46) 

t When the Jeans mass itself varies with time, linear gas fluctuations tend to be smoothed on 
a (filtering) scale that depends on the full thermal history of the gas instead of the instantaneous 
value of the sound speed (Gnedin & Hui 1998) . 
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where 

q ^ fe/Mpc-^ (3 47) 

This is the fit given by Bardeen et al. (1986) modified to account for the effects of baryon 
density following Sugiyama (1995). 

Another useful measurement of inhomogeneity is the quantity 

a^{M) = / dkk''P{k) \W{kR)\''. (3.48) 
271"^ Jo 

This is the mass variance of the density field smoothed by a window function W{kR) over 
a spherical volume of comoving radius R and average mass M = iJQf2mi?^/(2G). If the 
window is a top-hat in real space, then its Fourier transform is W{kR) = 3{kR)~^ (sin kR— 
kRcoskR). The significance of the window function is the following: the dominant con- 
tribution to cr(M) comes from perturbation components with wavelengths A = 27r/A; > R, 
because those with higher frequencies tend to be averaged out within the window vol- 
ume. As the fluctuation spectrum is falling with (lecrc;asing waves with much larger A 
contribute only a small amount. Hence, in terms of a mass M oc A^ oc k~^, we have 

cr{M) cx Af-(-''+")/". (3.49) 

Since P{k) is not a strict power law, n should be thought of as an approximate local value 
din P/ dink in the relevant range. For n > — 3, the variance <tm decreases with increasing 
M; this implies that, on the average, smaller masses condense out earlier than larger 
masses. Structures grow by the gradual separation and collapse of progressively larger 
units. Each parent unit will in general be made up of a number of smaller progenitor 
clumps that had collapsed earlier. This leads to a hierarchical pattern of clustering. In 
CDM, n > —3 at small M, increases with increasing M and reaches the asymptotic value 
n = 1 for M > 10-^^ M©. The power spectrum on galactic scales can be approximated 
by n w —2. 

3.3. Spherical collapse 

At late times the density contrast at a given wavelength becomes comparable to unity, 
linear perturbation theory fails at this wavelength, and a different model must be used 
to follow the collapse of bound dark matter systems ("halos"). Small scales are the first 
to become non-linear. Consider, at some initial time U, a spherical region of size r, and 
mass M that has a slight constant overdensity Si relative to the background pb, 

p{ti)=Pbiti)il + Si). (3.50) 

As the universe expands, the overdense region will expand more slowly compared to the 
background, will reach a maximum radius, and eventually collapse under its own gravity 
to form a bound virialized system. Such a simple model is called a spherical top-hat. As 
long as radial shells do not cross each other during the evolution, the motion of a test 
particle at (physical) radius r is governed by the equation (ignoring the vacuum energy 
component) 

d^r GM 

where M = (47r/3)r|/9b(<i)(l + Si) =const. Integrating we obtain 

1 /drV GM 



dt J 



=E, (3.52) 
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where £^ is a constant of integration. If > then will never become zero and the 
shell will expand forever. If < instead, then as r increases r will eventually become 
zero and later negative, implying a contraction and a collapse. Let's choose ti to be the 
time at which 5i is so small that the overdense region is expanding along with the Hubble 
flow. Then fi = {a/a)ri = H{ti)ri = HiVi at time ti, and the initial kinetic energy will 
be 

The potential energy at t = f j is 

\U\ = (^^) = G^p,{ti){l + (5,)r2 = ^Hfr^niil + <5,) = + 5,), (3.54) 

with n, dcnioting the initial value of the matter density parameter of the smooth back- 
ground universe. The total energy of the shell is therefore 

E = K,- K,n,il + S,) = K,n,inr^ - 1 - ,5,). (3.55) 

The condition < for the shell to eventually collapse becomes (1 + 6i) > or 

6^ > - 1. (3.56) 

In a Einstein-de Sitter universe at early times (f2j = 1), this condition is satisfied by any 
overdense region with Si > 0. In this case the patch will always collapse. 

Consider now a shell with i? < in a background Einstcin-de Sitter universe that 
expands to a maximum radius rmax ('turnaround', r = 0) and then collapses. The 
solution to the equation of motion can be written as the parametric equation for a 
cycloid 

r-max . , 6» - sin 6* 

r = — — (1 - COS^), i = tmax • (3.57) 

Z TT 

The background density evolves as pb{t) = {6TrGt'^)~^, and the density contrast becomes 

p{t) _ ^ _ 9 (6'-sin6 >)2 
Pbit) ^2(l-cos( 
For comparison, linear theory (in the limit of small t) gives 

This all agrees with what we knew already: at early times the sphere expands with the 
a DC f2/3 Hubble flow and density perturbations grow proportional to a. We can now see 
how linear theory breaks down as the perturbation evolves. There are three interesting 
epochs in the final stage of its development, which we can read directly from the above 
solutions. 



S=J^-l = - 1. (3.58) 



Turnaround. The sphere breaks away from the Hubble expansion and reaches a maxi- 
mum radius at 9 = n, t = tmax- At this point the true density enhancement with respect 
to the background is 

6=~-l = ^n'-l = 4.55, (3.60) 
which is definitely in the non-linear regime. By comparison, linear theory predicts 

= = 1-062. (3.61) 
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Collapse. If only gravity operates, then the sphere will collapse to a singularity at 
9 = 2^, t = 2tmax- This occurs when 

6l = ^(12t)'/' = 1.686. (3.62) 

Thus from equation (3.42) we see that a top-hat collapses at redshift z if its linear 
overdensity extrapolated to the present day is 

^c{z) = ^^=l.m^ + z), (3.63) 

where the second equality holds in a flat Einstcin-dc Sitter universe. This is termed the 
critical overdensity for collapse. An object of mass M collapsing at redshift z has an 
overdensity that is Vc times the linearly extrapolated density contrast today, ao{M), on 
that scale, 

i^c = 5c{z)/aa{M). (3.64) 

The mass scale associated with typical Vc = l-f non-linear fluctuations in the density 
field decreases from 10^^ Mq today to about 10^ at redshift 5. A.i z = 10, halos of 
10^° Mq collapse from much rarer = 3-cr fluctuations. 

Virialization. Collapse to a point at ^ = 27r (t = 2traa.^ = *vir) will never occur in 

practice as, before this happens, the approximation that matter is distributed in spherical 
shells and that the random velocities of the particles are small will break down. The dark 
matter will reach virial equilibrium by a process known as 'violent relaxation': since the 
gravitational potential is changing with time, individual particles do not follow orbits 
that conserve energy, and the net effect is to widen the range of energies available to 
them. Thus a time-varying potential can provide a relaxation mechanism that operates 
on the dynamical timescale rather than on the much longer two-body relaxation time. 
This process will convert the kinetic energy of collapse into random motions. 

At t = fmax all the energy is in the form of gravitational potential energy, and E = 
U = -GM/ryna^. At virialization U = -2K (virial theorem) and E = U + K = U/2 = 
— G'M/(2rvir). Hence r^ir — ?'inax/2. The mean density of the virialized object is then 
Pvir = 2^ Pmax, where Pmax is the density of the shell at turnaround. From equation 

(3.60) we have pmax = (9/16)7r^/9&(tmax), and P6(imax) = P6(tvir)(tvir/imax)^ = 4^pb{tvii)- 

Combining these relations we get: 

Pvir = 23p„^, = 23(9/16)^Vb(imax) = 2^ {9 / 16)TrH pi,{t^i,) = 18^2 (3 65) 
Therefore the density contrast at virialization in an Einstein-de Sitter universe is 

Ac = 178. (3.66) 

In a universe with a cosmological constant, the collapse of a top-hat spherical perturba- 
tion is described by 

GM A 

^=-— +r- ^'-''^ 

In a flat + ^a = I cosmology, the final overdensity relative to the critical density gets 
modified according to the fitting formula (Bryan & Norman 1998) 

Ac = 18 7r^-|-82c^-39c^^ (3.68) 

where d = nf„ — 1 is evaluated at the collapse redshift. A spherical top- hat collapsing 
today in a universe with Qm = 0.24 has a density contrast at virialization of Ac = 
93. This corresponds to an overdensity relative to the background matter density (= 
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i^mPc) of Ac/r^TO = 388: the faster expansion of a low-density universe means that 
the perturbation turns around and collapses when a larger density contrast has been 
produced. For practical reason a density contrast of 200 relative to the background is 
often used to define the radius, r2005 that marks the boundary of a virialized region. 

A halo of mass M collapsing at redshift z can be described in terms of its virial radius 
Tvir, circular velocity Vc and virial temperature Tvir (Barkana & Loeb 2001): 



2GM 



._-g2..e..o^K(^)''%v,(i±i), ,,n, 

where / = {Q,m/^m){^c/^8''^'^) and fj. = 0.59 is the mean molecular weight for fully 
ionized primordial gas (/z = 1.23 for neutral primordial gas). The binding energy of the 
halo is approximately 

where the coefRcient of 1/2 is exact for a singular isothermal sphere. The binding energy 
of the baryons is smaller by a factor Q^/^b — 5.7. The energy deposition by supcrnovae 
in the shallow potential wells of subgalactic systems may then lift out metal-enriched 
material from the host (dwarf) halos, causing the pollution of the IGM at early times 
(Madau et al. 2001). 



3.4. Dark halo mergers 
The assumption that virialized objects form from smooth spherical collapse, while pro- 
viding a useful framework for thinking about the formation histories of gravitationally 
bound dark matter halos, does not capture the real nature of structure formation in 
CDM theories. In these models galaxies are assembled hierarchically through the merg- 
ing of many smaller subunits that formed in a similar manner at higher redshift (see Fig. 
5). Galaxy halos experience multiple mergers during their lifetime, with those between 
comparable-mass systems ("major mergers") expected to result in the formation of el- 
liptical galaxies (see, e.g., Barnes 1988; Hernquist 1992). Figure 6 shows the number 
of major mergers (defined as mergers where the mass ratio of the progenitors is > 0.3) 
per unit redshift bin experienced by halos of different masses. For galaxy-sized halo this 
quantity happens to peak in the redshift range 2-4, the epoch when the observed space 
density of optically-selected quasar also reaches a maximum. 

The merger between a large parent halo and a smaller satellite system will evolve 
under two dynamical processes: dynamical friction, that causes the orbit of the satellite 
to decay toward the central regions, and tidal stripping, that removes material from the 
satellite and adds it to the diffuse mass of the parent. Since clustering is hierarchical, 
the satellite will typically form at earlier times and have a higher characteristic density 
and a smaller characteristic radius. The study of the assembly history of dark matter 
halos by repeated mergers is particularly well suited to the N-body methods that have 
been developed in the past two decades. Numerical simulations of structure formation by 
dissipationless hierarchical clustering from Gaussian initial conditions indicate a roughly 
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Figure 5. An example of a merger tree obtained from an N-body simulation of a 9 x W^^h^^ Mq 
halo at redsliift z = 0. Each circle represents a dark matter halo identified in the sinmlation, 
the area of the circle being proportional to halo mass. The vertical position of each halo on the 
plot is determined by the redshift z at which the halo is identified, the horizontal positioning 
is arbitrary. The solid lines connect halos to their progenitors. The solid line in the panel on 
the left-hand side shows the fraction of the final mass contained in resolved progenitors as a 
function of redshift. The dotted line shows the fraction of the final mass contained in the largest 
progenitor as a function of redshift. (Prom Helly et al. 2003.) 



universal spherically-averaged density profile for the resulting halos (Navarro et al. 1997, 
hereafter NFW): 

pNFwir) = — — —, (3.73) 

ca;(l -I- ca;)^ 

where x = r/r^i-c and the characteristic density ps is related to the concentration param- 
eter c by 

SttG 3 ln(l + c)-c/(l + c)' ^ ' 

This function fits the numerical data presented by NFW over a radius range of about two 
orders of magnitude. Equally good fits arc obtained for high-mass (rich galaxy cluster) 
and low-mass (dwarf) halos. Power-law fits to this profile over a restricted radial range 
have slopes that steepen from —1 near the halo center to —3 at large cr/rvir. Bullock 
et al. (2001) found that the concentration parameter follows a log-normal distribution 
where the median depends on the halo mass and redshift, 

c^ed(M,^) = -^f^V, (3.75) 



1 + z VM*_ 

where M* is the mass of a typical halo collapsing today. The halos in the simulations by 
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Figure 6. Mean number of major mergers experienced per unit redshift by halos with masses 
> 10^° M0. Solid line: progenitors of a Mo = 2 x 10^^ Mq halo at « = 0. Dashed line: same 
for Mo = 3 X 10^^ M0. Dotted line: same for Mq = lO" Mq. (Prom Volonteri et al. 2003.) 

Bullock et al. were best described by a = —0.13 and = 9.0, with a scatter around the 
median of a\ogc = 0.14 dex. 

"Via Lactea" , the highest resolution N-body simulation to date of the formation of a 
Milky Way-sized halo (Diemand et al. 2007a, b), shows that the fitting formula proposed 
by NFW with concentration c = 12 provides a reasonable approximation to the density 
profile down a convergence radius of rconv = 1-3 kpc. Within the region of convergence, 
deviations from the best-fit NFW matter density are typically less than 10%. From 10 
kpc down to rconv Via Lactea is actually denser than predicted by the NFW formula. 
Near rconv the density approaches the NFW value again while the logarithmic slope is 
shallower (—1.0 at rconv) than predicted by the NFW fit. 

3.5. Assemhly history of a Milky Way halo 

The simple spherical top-hat collapse ignores shell crossing and mixing, accretion of self- 
bound clumps, triaxiality, angular momentum, random velocities, and large scale tidal 
forces. It is interesting at this stage to use "Via Lactea" and study in more details, 
starting from realistic initial conditions, the formation history of a Milky Way-sized halo 
in a ACDM cosmology. The Via Lactea simulation was performed with the PKDGRAV 
tree-code (Stadel 2001) using the best-fit cosmological parameters from the WMAP 3- 
year data release. The galaxy-forming region was sampled with 234 million particles of 
mass 2.1 x 10^ M©, evolved from redshift 49 to the present with a force resolution of 90 
pc and adaptive time-steps as short as 68,500 yr, and centered on an isolated halo that 
had no major merger after z = 1.7, making it a suitable host for a Milky Way-like disk 
galaxy (see Fig. 8). The number of particles is an order of magnitude larger than used in 
previous simulations. The rmi was completed in 320,000 CPU hours on NASA's Project 
Columbia supercomputer, currently one of the fastest machines available. (More details 
about the Via Lactea run are given in Diemand et al. 2007a, b. Movies, images, and data 
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Figure 7. Top: logarithmic slope of the density profile of the "Via Lactea" halo, as a function 
of radius. Densities were computed in 50 radial logarithmic bins, and the local slope was 
determined by a finite difference approximation using one neighboring bin on each side. The 
thin line shows the slope of the best-fit NFW profile with concentration c = 12. The vertical 
dotted line indicates the estimated convergence radius: local densities (but not necessarily the 
logarithmic slopes) should be correct to within 10% outside of this radius. Bottom: the residuals 
in percent between the density profile and the best-fit NFW profile, as a function of radius. Here 
'"200 is the radius within which the enclosed average density is 200 times the background value, 
'"200 ~ 1.35 rvir in the adopted cosmology. (From Diemand et al. 2007a.) 



are available at http://www.ucolick.org/'^diemand/vl.) The host halo mass at z = is 
M200 = 1.8 X 10^^ M0 within a radius of r2oo = 389 kpc. 

Following the spherical top-hat model, the common procedure used to describe the 
assembly of a dark matter halo is to define at each epoch a virial radius rvir (or, equiv- 
alently, r2oo), which depends on the cosmic background density at the time. As the 
latter decreases with the Hubble expansion, formal virial radii and masses grow with 
cosmic time even for stationary halos. Studying the transformation of halo properties 
within Tvir (or some fraction of it) mixes real physical change with apparent evolutionary 
effects caused by the growing radial window, and makes it hard to disentangle between 
the two. Figure 9 shows the formation of Via Lactea where radial shells enclosing a 
fixed mass, Tm, have been used instead. Unlike Tvir, ?"m stops growing as soon as the 
mass distribution of the host halo becomes stationary on the corresponding scale. Note 
that mass and substructure are constantly exchanged between these shells, as tm is not a 
Lagrangian radius enclosing the same material at all times, just the same amount of it. 
The fraction of material belonging to a given shell in the past that still remains within 
the same shell today is shown in Figure 10. The mixing is larger before stabilization, 
presumably because of shell crossing during collapse, and smaller near the center, where 
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Figure 8. Projected dark matter density-squared map of our simulated Milky Way halo region 
("Via Lactea") at different redshift, from z = 10.6 to the present. The image covers an area of 
400 X 300 physical kpc. 



most of the mass is in a dynamically cold, concentrated old component (Diemand et al. 
2005b). Outer shells number 9 and 10, for example, retain today less than 25% of the 
particles that originally belonged to them at a < 0.4. Note that the collapse times also 
appear to differ from the expectations of spherical top-hat. Shell number five, for ex- 
ample, encloses a mean density of about 100 pj! today, a virial mass of 1.5 x 10^^ Mq 
and should have virialized just now according to spherical top-hat. It did so instead 
much earlier, at a = 0.6. Even the next larger shell with 1.8 x 10^^ Mq stabilized before 
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Figure 9. Evolution of radii tm enclosing a fixed mass versus cosmic time or scale factor a. 
The enclosed mass grows in constant amounts of 0.3 x 10^^ M© from bottom to top. Initially all 
spheres are growing in the physical (uoii comoving) units used here. Inner shells turn around, 
collapse and stabilize, while the outermost shells are still expanding today. Solid circles: points 
of maximum expansion at the turnaround time tmax- Open squares: time after turnaround 
where rm first contracts within 20% of the final value. These mark the approximate epoch of 
"stabilization" . 



a = 0.8. It appears that spherical top-hat provides only a crude approximation to the 
virialized regions of simulated galaxy halos. 

To understand the mass accretion history of the Via Lactea halo it is useful to analyze 
the evolution of mass within fixed physical radii. Figure 11 shows that the mass within 
all radii from the resolution limit of ~1 kpc up to 100 kpc grows during a series of major 
mergers before a = 0.4. After this phase of active merging and growth by accretion the 
halo mass distribution remains almost perfectly stationary at all radii. Only the outer 
regions (~ 400 kpc) experience a small amount of net mass accretion after the last major 
merger. The mass within 400 kpc increases only mildly, by a factor of 1.2 from z = 1 to 
the present. During the same time the mass within radii of 100 kpc and smaller, the peak 
circular velocity and the radius where this is reached, all remain constant to within 10%. 
The fact that mass definitions inspired by spherical top-hat fail to accurately describe 
the real assembly of galaxy halos is clearly seen in Figure 11, where M200 is shown to 
increase at late times even when the halo physical mass remains the same. This is just 
an artificial effect caused by the growing radial windows rvir and r2oo as the background 
density decreases. For Via Lactea M200 increases by a factor of 1.8 from z = 1 to the 
present, while the real physical mass within a 400 kpc sphere grows by only a factor of 
1.2 during the same time interval, and by an even smaller factor at smaller radii. 

3.6. Smallest SUSY- CDM microhalos 
As already mentioned above, the key idea of the standard cosmological paradigm for 
the formation of structure in the universe, that primordial density fluctuations grow 
by gravitational instability driven by cold, coUisionless dark matter, is constantly being 
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Figure 10. Fraction of material belonging to shell i at epoch a that remains in the same shell 
today. Shells are same as in Fig. 9, numbered from one (inner) to ten (outer). Solid circles: time 
of maximum expansion. Open squares: stabilization epoch. Mass mixing generally decreases 
with time and towards the halo center. 

elaborated upon and explored in detail through supercomputer simulations, and tested 
against a variety of astrophysical observations. The leading candidate for dark matter 
is the neutralino, a weakly interacting massive particle predicted by the supersymmetry 
(SUSY) theory of particle physics. While in a SUSY-CDM scenario the mass of bound 
dark matter halos may span about twenty order of magnitudes, from the most massive 
galaxy clusters down to Earth-mass clumps (Green et al. 2004) , it is the smallest micro- 
halos that collapse first, and only these smallest scales are affected by the nature of the 
relic dark matter candidate. 

Recent numerical simulations of the collapse of the earliest and smallest gravitationally 
bound CDM clumps (Dicmand et al. 2005a; Gao et al. 2005) have shown that tiny 
virialized microhalos form at redshifts above 50 with internal density profiles that are 
quite similar to those of present-day galaxy clusters. At these epochs a significant fraction 
of neutralinos has already been assembled into non-linear Earth-mass overdensities. If 
this first generation of dark objects were to survive gravitational disruption during the 
early hierarchical merger and accretion process - as well as late tidal disruption from 
stellar encounters (Zhao et al. 2007) - then over 10"'^^ such clumps may populate the halo 
of the Milky Way. The nearest microhalos may be among the brightest sources of 7-rays 
from neutralino annihilation. As the annihilation rate increases quadratically with the 
matter density, small-scale dumpiness may enhance the total 7-ray flux from nearby 
extragalactic systems (like M31), making them detectable by the forthcoming GLAST 
satellite or the next-generation of air Cerenkov telescopes. 

The possibility of observing the fingerprints of the smallest-scale structure of CDM in 
direct and indirect dark matter searches hinges on the ability of microhalos to survive 
the hierarchical clustering process as substructure within the larger halos that form at 
later times. In recent years high-resolution N-body simulations have enabled the study 
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a = 1 / (1+z) 



Figure 11. Mass accretion history of Via Lactea. Masses within spheres of fixed physical radii 
centered on the main progenitor are plotted against the cosmological expansion factor a. The 
thick solid lines correspond to spheres with radii given by the labels on the left. The thin solid 
lines correspond to nine spheres of intermediate radii that are 1.3, 1.6, 2.0, 2.5, 3.2, 4.0, 5.0, 
6.3 and 7.9 times larger than the next smaller labeled radius. Dashed line: M2oo- The halo is 
assembled during a phase of active merging before a ~ 0.37 {z ~ 1.7) and its net mass content 
remains practically stationary at later times. 



of gravitationally-bound subhalos with Msuh/M > 10~^ on galaxy (and galaxy cluster) 
scales (e.g. Moore et al. 1999; Klypin et al. 1999; Stoehr et al. 2003). The main differences 
between these subhalos - the surviving cores of objects which fell together during the 
hierarchical assembly of galaxy-size systems - and the tiny sub-microhalos discussed here 
is that on the smallest CDM scale the effective index of the linear power spectrum of 
mass density fluctuations is close to —3. In this regime typical halo formation times 
depend only weakly on halo mass, the capture of small clumps by larger ones is very 
rapid, and sub-microhalos may be more easily disrupted. 

In Diemand et al. (2006) we presented a large N-body simulation of early substructure 
in a SUSY-CDM scenario characterized by an exponential cutoff in the power spectrum 
at lO"*^ Mq. The simulation resolves a 0.014 Mq parent "SUSY" halo at z = 75 with 
14 million particles. Compared to a z = galaxy cluster substructure within the SUSY 
host is less evident both in phase-space and in physical space (see Figs. 12 and 13), and 
it is less resistant against tidal disruption. As the universe expands by a factor of 1.3, 
between 20 and 40 percent of well-resolved SUSY substructure is destroyed, compared 
to only ~ 1 percent in the low-redshift cluster. Nevertheless SUSY substructure is just 
as abundant as in z = galaxy clusters, i.e. the normalized mass and circular velocity 
functions are very similar. 
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Figure 12. Local dark matter density maps. The left panels illustrate the almost simultaneous 
structure formation in the SUSY run at different epochs, within a sphere of physical radius 
r including a mass of 0.014 Mq. A galaxy cluster halo {right panels) forms in the standard 
hierarchical fashion: the dark matter distribution is shown within a sphere of radius r including 
a mass of 5.9 x 10^* Mq. The SUSY and cluster halos have concentration parameter for a NFW 
profile of c = 3.7 and c = 3.5, respectively. In each image the logarithmic color scale ranges 
from 10 to 10^ pc{z). 

4. The Dawn of Galaxies 

4.1. Uncertainties in the power spectrum 

As mentioned in the Introduction, some shortcomings on galactic and sub-galactic scales 
of the currently favored model of hierarchical galaxy formation in a universe dominated 
by CDM have recently appeared. The significance of these discrepancies is still debated, 
and 'gastrophysical' solutions involving feedback mechanisms may offer a possible way 
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Figure 13. Phase-space density (p/a^, where a is the one-dimensional velocity dispersion) map 
for a z = 75 SUSY halo (left) and a z = galaxy cluster (right). Note the diflerent color scales: 
relative to the average phase-space density, the logarithmic color scale ranges from 10 to 10^ in 
the SUSY halo and from 10 to 10^ in the cluster halo. 

out. Other models have attempted to solve the apparent small-scale problems of CDM 
at a more fundamental level, i.e. by reducing small-scale power. Although the 'standard' 
ACDM model for structure formation assumes a scale-invariant initial power spectrum 
of density fluctuations with index n = 1, the recent WMAP data find strong evidence for 
a departure from scale invariance, with a best fit value n = 0.951^Q Qi[g. Furthermore the 
1-year WMAP results favored a slowly varying spectral index, dn/d\nk — — O.OSllQpl^g, 
i.e. a model in which the spectral index varies as a function of wavenumber k (Spergel 
et al. 2003). In the 3-year WMAP data such a 'running spectral index' leads to a marginal 
improvement in the fit. Models with either n < 1 or dn/dlnfc < predict a significantly 
lower amplitude of fluctuations on small scales than standard ACDM. The suppression of 
small-scale power has the advantage of reducing the amount of substructure in galactic 
halos and makes small halos form later (when the universe was less dense) hence less 
concentrated (Zentner & Bullock 2002). 

Figure 14 shows the linearly extrapolated (to z = 0) variance of the mass-density field 
for a range of cosmological parameters. Note that the new WMAP results prefer a low 
value for erg, the rms mass fluctuation in a 8 ft." ^ Mpc sphere. This is consistent with a 
normalization by the z = X-ray cluster abundance (Reiprich & Boringer 2002). For 
comparison we have also included a model with a higher normalization, ag = 0.9, from 
the best-fit model to 1-year WMAP data. As in the CDM paradigm structure formation 
proceeds bottom-up, it then follows that the loss of small-scale power modifies structure 
formation most severely at the highest redshifts, significantly reducing the number of 
self-gravitating objects then. This, of course, will make it more difficult to reionize the 
universe early enough. 

It has been argued, for example, that one popular modification of the CDM paradigm, 
warm dark matter (WDM), has so little structure at high redshift that it is unable to 
explain the WMAP observations of an early epoch of reionization (Barkana et al. 2001). 
And yet the WMAP running-index model may suffer from a similar problem. A look 
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Figure 14, The variance of the matter-density field vs. mass M, for several different cosmolo- 
gies, all based on WMAP results. Solid curve: 3-year WMAP-only best-fit model. Dotted curve: 
AWDM with a particle mass mx ~ 2 keV, otherwise same as before. Dashed curve: 3-year year 
WMAP best-fit running spectral index model. Dash-dotted curve: 1-year WMAP-only best-fit 
tilted model. Here n refers to the spectral index at fc = 0.05 Mpc~^. The horizontal line at the 
top of the figure shows the value of the extrapolated critical collapse overdensity Sc{z) at the 
reionization redshift z = 11. 



at Figure 14 shows that 10^ M0 halos will collapse at z = 11 from 2,4 ct fluctuations in 
a tilted ACDM model with n = 0.951 and as = 0.74 (best-fit 3-year WMAP model), 
but from much rarer 3.7 and 3.6 a fluctuations in the WDM and running-index model, 
respectively. The problem is that scenarios with increasingly rarer halos at early times 
require even more extreme assumptions (i.e. higher star formation efficiencies and UV 
photon production rates) in order to be able to reionize the universe suitably early (e.g. 
Somerville et al, 2003; Wyithe & Loeb 2003; Ciardi et al, 2003; Cen 2003), Figure 15 de- 
picts the mass fraction in all collapsed halos with masses above the cosmological filtering 
mass for a case without reionization and one with reionization occurring at 2; ~ 11, At 
early epochs this quantity appears to vary by orders of magnitude in different models, 

4,2, First baryonic objects 

The study of the non-linear regime for the baryons is far more complicated than that of 
the dark matter because of the need to take into account pressure gradients and radiative 
processes. As a dark matter halo grows and virializes above the cosmological Jeans mass 
through merging and accretion, baryonic material will be shock heated to the effective 
virial temperature of the host and compressed to the same fractional overdensity as the 
dark matter. The subsequent behavior of gas in a dark matter halo depends on the 
efficiency with which it can cool. It is useful here to identify two mass scales for the host 
halos: (1) a molecular cooling mass Mn^ above which gas can cool via roto-vibrational 
levels of H2 and contract, M^^ ~ 10^ [(1 -f z)/10]~^/^ Mq (virial temperature above 
200 K); and (2) an atomic cooling mass Afjj above which gas can cool efficiently and 
fragment via excitation of hydrogen Lya, Mh ~ 10® [(1 -I- z)/10]~^/^ Mq (virial temper- 
ature above 10** K). Figure 16 shows the cooling mechanisms at various temperatures for 
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Figure 15. Mass fraction in all collapsed halos above the cosmological filtering (Jeans) mass 
as a function of redshift, for different power spectra. Curves are the same as in Figure 14. 
Left panel: filtering mass computed in the absence of reionization. Right panel: filtering mass 
computed assuming the universe is reionized and reheated to Te — 10* K by ultraviolet radiation 
at 2: ~ 11. 
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Figure 16. Cooling rates for Bremsstrahlung (dotted line), H (dashed line) and He (dash- dotted 
line) line cooling, and H2 (circles) cooling. The e~ , H 11, He 11, and He iii abundances were com- 
puted assuming colhsional equilibrium, and the H2 fractional abundance was fixed at 3 x 10~*, 
which is typical for early objects. At temperatures between 100 and 10,000 K, the H2 molecule 
is the most effective coolant (From Fuller & Couchman 2000). 



primordial gas, while Figure 17 shows the fraction of the total mass in the universe that 
is in collapsed dark matter halos with masses greater than M-^^ and il/n at different 
epochs. 

High-resolution hydrodynamics simulations of early structure formation are a powerful 
tool to track in detail the thermal and ionization history of a clumpy IGM and guide 
studies of primordial star formation and reheating. Such simulations performed in the 
context of ACDM cosmologies have shown that the first stars (the so-called 'Population 
III') in the universe formed out of metal-free gas in dark matter minihalos of mass above 
a few X 10^ M0 (Abel et al. 2000; Fuller & Couchman 2000; Yoshida et al. 2003; Reed 
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Figure 17. Sohd lines: Total mass fraction in all collapsed dark matter halos above the 
molecular cooling and the atomic cooling masses, Mb_2 and Mh, as a function of redshift. 

et al. 2005) condensing from the rare high-j^c peaks of the primordial density fluctuation 
field at z > 20, and were likely very massive (e.g. Abel et al. 2002; Bromm et al. 2002; 
see Bromm & Larson 2004 and Ciardi & Ferrara 2005 for recent reviews). In Kuhlen 
& Madau (2005) we used a modified version of enzo, an adaptive mesh refinement 
(AMR), grid-based hybrid (hydro+N-body) code developed by Bryan & Norman (see 
http://cosmos.ucsd.edu/enzo/) to solve the cosmological hydrodynamics equations and 
study the cooling and collapse of primordial gas in the first baryonic structures. The 
simulation samples the dark matter density field in a 0.5 Mpc box with a mass resolution 
of 2000 Mq to ensure that halos above the cosmological Jeans mass are well resolved at 
all redshifts z < 20. The AMR technique allows us to home in, with progressively finer 
resolution, on the densest parts of the "cosmic web" . During the evolution from z = 99 to 
z — 15, refined (child) grids are introduced with twice the spatial resolution of the coarser 
(parent) grid when a cell reaches a dark matter overdensity (baryonic overdensity) of 2.0 
(4.0). Dense regions are allowed to dynamically refine to a maximum resolution of 30 
pc (comoving). The code evolves the non-equilibrium rate equations for nine species (H, 
H+, H^, e. He, He"'", He"*""*", H2, and H2''") in primordial self-gravitating gas, including 
radiative losses from atomic and molecular line cooling, and Compton cooling by the 
cosmic background radiation. 

The clustered structure around the most massive peaks of the density field is clearly 
seen in Figure 18, a 3D volume rendering of the simulated volume at redshifts 21 and 15.5. 
The figure shows gas at 4 < (5b < 10, with the locations of dark matter minihalos marked 
by spheres colored and sized according to their mass (the spheres are only markers, the 
actual shape of the halos is typically non-spherical). Several interleaving filaments are 
visible, at the intersections of which minihalos are typically found. At z = 21, 55 bound 
halos are indentified in the simulated volume: by z = 17.5 this number has grown to 149, 
and by z = 15.5 to 262 halos. At this epoch, only four halos have reached the critical 
virial temperature for atomic cooling, Tvir = lO'' K. 
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Figure 18. A 3D volume rendering of the IGM in the inner 0.5 Mpc simulated box at 2 = 21 
(lejt panetj and z — 15.5 {right panel). Only gas with overdensity 4 < 5;, < 10 is shown: the 
locations of dark matter minihalos are marked by spheres with sizes proportional to halo mass. 
At these epochs, the halo finder algorithm identifies 55 {z = 21) and 262 {z = 15.5) bound 
clumps within the simulated volume. (From Kuhlen & Madau 2005.) 
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Figure 19. Lejt: Fraction of cold and cold+dense gas within the virial radius of all halos 
indentified at 2 = 17.5 with Tvir > 400 K, as a function of halo mass. Triangles: fc, fraction of 
halo gas with T < 0.5 Tvir and overdensity Si, > 1000 that cooled via roto- vibrational transitions 
of H2. Diamonds: fed, fraction of gas with T < O.STVir and p > 10^^ M0 IVIpc"'^ that is 
available for star formation. The straight lines represent mean regression analyses of fc and 
/cd with the logarithm of halo mass. Right: Mass-weighted mean H2 fraction as a function of 
virial temperature for all halos at 2 = 17.5 with Tvir > 400 K and f^d < 0.1 (empty circles) 
or fed > 0.1 {filled circles). The straight line marks the scaling of the temperature-dependent 
asymptotic molecular fraction. 



The primordial fractional abundance of H2 in the low-density IGM is small, xb_2 = 
[H2/H] ~ 2 X 10~^, as at 2; > 100 H2 formation is inhibited because the required inter- 
mediaries, either or H^, are destroyed by CIvIB photons. IVIost of the gas in the 
simulation therefore cools by adiabatic expansion. Within collapsing minihalos, however, 
gas densities and temperatures are large enough that H2 formation is catalyzed by H~ 
ions through the associative detachment reaction H-|-H~ — > H2-fe^, and the molecular 
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fraction increases at the rate dxa^/dt oc XenmT^ir^ , where Xe is the number of electrons 
per hydrogen atom. For Tvir ^ a few thousand kelvins the virialization shock is not ion- 
izing, the free electrons left over from recombination are depleted in the denser regions, 
and the production of H2 stalls at a tcmpcratmo-dcpcndcnt asymptotic molecular frac- 
tion ~ ^^~^Tlir I'^Cl + t/trcc), where tree IS thc hydrogcn recombination time-scale 
(Tegmark et al. 1997). A typical H2 fraction in excess of 200 times thc primordial value 
is therefore produced after the collapse of structures with virial temperatures of order 
10^ K. This is large enough to efficiently cool the gas and allow it to collapse within a 
Hubble time unless significant heating occurs during this phase (Abel et al. 2000; Yoshida 
et al. 2003). 

Fig. 19 (left panel) shows the fraction of cold gas within the virial radius as a function 
of halo mass for all the halos identified at redshift 17.5. Following Machacck ct al. 
(2001), we define /c as the fraction of gas with temperature < O.STvir and density 
> 1000 times the background (this is the halo gas that is able to cool below the virial 
temperature because of H2), and fed as the fraction of gas with temperature < 0.5 Tvb 
and (physical) density > 10^^ Mq Mpc~^ (this is the self-gravitating gas available for star 
formation). As in Machacek et al. (2001), we find that both fc and fed are correlated 
with halo mass. The threshold for significant baryonic condensation (non-zero /cd) is 
approximately 5 x 10^ Mq at these redshifts (Haiman et al. 1996). Also depicted in 
Fig. 19 (right panel) is the mass-weighted mean molecular fraction of all halos with 
Tvir > 400 K. Filled circles represent halos with fed > 0.1, while open circles represent 
the others. The straight line marks the scaling of the asymptotic molecular fraction in 
the electron-depletion transition regime. The maximum gas density reached at redshift 
15 in the most refined region of our simulation is 4 x 10"'' cm^'^ (corresponding to an 
over density of 3 x 10*): within this cold pocket the excited states of H2 are in LTE and 
the cooling time is nearly independent of density. 

4.3. 21cm signatures of the neutral IGM 

It has long been known that neutral hydrogen in the diffuse IGM and in gravitationally 
collapsed structures may be directly detectable at frequencies corresponding to thc red- 
shifted 21cm line of hydrogen (Field 1959; Sunyaev & Zel'dovich 1975; Hogan & Rees 
1979, see Furlanetto et al. 2006 for a recent review). The emission or absorption of 21cm 
photons from neutral gas is governed by the spin temperature T5, defined as 

m/no = 3exp(-T*/Ts). (4.76) 

Here no and ui arc thc number densities of hydrogen atoms in the singlet and triplet 
n = 1 hyperfine levels, and T* = 0.068 K is the temperature corresponding to the energy 
difference between the levels. To produce an absorption or emission signature against 
the CMB, thc spin temperature must differ from thc temperature; of the CMB, T = 
2.73 (1+2) K. At 30 < z < 200, prior to the appearance of non-linear baryonic objects, the 
IGM cools adiabatically faster than the CMB, spin-exchange collisions between hydrogen 
atoms couple T5 to thc kinetic temperature Tg of the cold gas, and cosmic hydrogen can 
be observed in absorption (Scott & Rees 1990; Loeb & Zaldarriaga 2004). At lower 
redshifts, the Hubble expansion rarefies the gas and makes collisions inefficient: the spin 
states go into equilibrium with thc radiation, and as T5 approaches T the 21cm signal 
diminishes. It is the first luminous sources that make uncoUapsed gas in the universe 
shine again in 21cm, by mixing the spin states either via Lya scattering or via enhanced 
free electron-atom collisions (e.g. Madau ct al. 1997; Tozzi ct al. 2000; Nusser 2005). 

While the atomic physics of the 21cm transition is well understood in the cosmological 
context, exact calculations of the radio signal expected during the era between the col- 
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lapse of the first baryonic structures and the epoch of complete reionization have been 
difficult to obtain, as this depends on the spin temperature, gas overdensity, hydrogen 
neutral fraction, and line-of-sight peculiar velocities. When Ts = Te, the visibility of the 
IGM at 21cm revolves around the quantity (Tg - T)/Te. If < T, the IGM will appear 
in absorption against the CMB; in the opposite case it will appear in emission. To deter- 
mine the kinetic temperature of the IGM during the formation of the first sources, one 
needs a careful treatment of the relevant heating mechanisms such as photoionization 
and shock heating. In addition to the signal produced by the cosmic web, minihalos with 
virial temperatures of a few thousand kelvins form in abundance at high redshift, and 
are sufficiently hot and dense to emit coUisionally-excited 21 cm radiation (Iliev et al. 
2002). 

The H I spin temperature is a weighted mean between Tg and T, 

Ts = ^^^^. (4.77) 

1 + 2/ 

where coupling efficiency y is the sum of three terms, 

y = ^(CH + Ce + Cp). (4.78) 

Here A = 2.85 x 10~^^s~^ is the spontaneous emission rate and Ch, Cg, and Cp are 
the de-excitation rates of the triplet due to collisions with neutral atoms, electrons, 
and protons. A fourth term must be added in the presence of ambient Lya radiation, 
as intermediate transitions to the 2p level can mix the spin states and couple Ts to 
Te, the "Wouthuysen-Field" effect (Wouthuysen 1952; Field 1958; Hirata 2006). In 
the absence of Lya photons, to unlock the spin temperature of a neutral medium with 
(say) Te = 500 K from the CMB requires a collision rate Ch > AT/T^, corresponding 
to a baryon overdensity St > 5 [(1 + z)/20]~'^. Not only dense gas within virialized 
minihalos but also intergalactic filamentary structure heated by adiabatic compression 
or shock heating may then be observable in 21cm emission. A population of X-ray sources 
( "miniquasars" ) turning on at early stages (e.g. Madau et al. 2004; Ricotti et al. 2005) 
may make even the low-density IGM visible in 21cm emission as structures develop in 
the pre-reionization era. 

In Kuhlen et al. (2006) we used the same hydrodynamical simulations of early structure 
formation in a ACDM universe discussed in § 1.4.2 to investigate the spin temperature 
and 21cm brightness of the diffuse IGM prior to the epoch of cosmic reionization, at 
10 < 2; < 20. The two-dimensional distribution of gas overdensity and spin temperature 
at z = 17.5 is shown in Figure 20 for a simulation with no radiation sources ("NS") and 
one ("MQ") in which a miniquasar powered by a 150 Mq black hole turns on at redshift 
21 within a host halo of mass 2 x 10^ M©. The miniquasar shines at the Eddington 
rate and emits X-ray radiation with a power- law energy spectrum oc E^^ in the range 
0.2—10 keV. The color coding in this phase diagram indicates the fraction of the simulated 
volume at a given {6b, Ts). In both runs we have assumed no Lya mixing, so that the 
visibility of hydrogen at 21cm is entirely determined by collisions. Only gas with neutral 
fraction > 90% is shown in the figure. The low-density IGM in the NS run lies on the 
yellow Ts = T = 50.4 K line: this is gas cooled by the Hubble expansion to <C T 
that cannot unlock its spin states from the CMB and therefore remains invisible. At 
overdensities between a few and ~ 200, H-H collisions become efficient and adiabatic 
compression and shocks from structure formation heats up the medium well above the 
radiation temperature. The coupling coefficient at this epoch is y ^ SbT~'^-^'^: gas in 
this regime has T < Tg ~ yTe < Te and appears in emission against the CMB (red and 
green swath). Some residual hydrogen with overdensity up to a few tens, however, is 
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Figure 20. Two-dimensional distribution of spin temperature versus baryonic overdensity at 
z = 17.5. The color coding indicates the fraction of the simulated volume at a given {5b, Ts). 
Left: NS run. The volume and mass-averaged spin temperatures are 48.6 and 67.5 K, respec- 
tively. Right: MQ run. Only gas with neutral fraction > 90% is shown in the figure. The volume 
and mass-averaged spin temperatures are 82.6 and 138.0 K, respectively. The red horizontal 
line marks the temperature of the CMB at that redshift. 

still colder than the CMB, and is detectable in absorption. At higher densities, y ^ 1 
and Tg Tf.: the blue cooling branch follows the evolutionary tracks in the kinetic 
temperature-density plane of gas shock heated to virial values, = 2000— 10"* K, which 
is subsequently cooling down to 100 K because of H2 line emission. 

The effect of the miniquasar on the spin temperature is clearly seen in the right panel 
of Figure 20. X-ray radiation drives the volume-averaged temperature and electron 
fraction (xg) within the simulation box from (8 K, 1.4 x 10"'*) to (2800 K, 0.03), therefore 
producing a warm, weakly ionized medium (Kuhlen & Madau 2005). The H-H collision 
term for spin exchange in the low-density IGM increases on the average by a factor 
350" '^^ ~ 8, while the e-H collision term grows to Ce ~ 0.5Ch. Gas with {Si,,Te,Xe) = 
(1, 2800 K, 0.03) has coupling efhciency y = 0.008, spin temperature = 73K > T, 
and can now be detected in emission against the CMB. Within 150 comoving kpc from 
the source, the volume-averaged electron fraction rises above 10%, and e-H collisions 
dominates the coupling. 

A beam of 21cm radiation passing through a neutral hydrogen patch having opti- 
cal depth T and spin temperature T5 causes absorption and induces emission. In the 
comoving frame of the patch, the radiative transfer equation yields for the brightness 
temperature through the region: Tf, — Ter'^ + Ts{l — e~'^). We have used our numer- 
ical simulations to perform 21cm radiation transport calculations, including the effect 
of peculiar velocities and local changes in spin temperature, gas density, and neutral 
hydrogen fraction. The resulting 21cm radio signal is shown in Figure 21 (lower panels), 
together with an image of the projected hydrogen spin temperature (upper panels): the 
latter highlights the abundance of structure within our simulation box on scales up to 
hundreds of kpc. Due to Hubble and peculiar velocity shifts not all of this structure 
contributes to the STi, map. In the NS simulation, coherent features in the IGM can be 
discerned in emission (Ts > T, ST), > 0) : this filamentary shock-heated structure is typ- 
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Figure 21. Projected (mass- weighted) spin temperature {upper panels, logarithmic scale) and 
21cm differential brightness temperature {lower panels, composite linear scale) in a 0.5 Mpc 
simulation box for runs "NS" {left) and "MQ" {right) at z = 17.5. The location of the miniquasar 
is indicated by crosses in the right panels. 



ically surrounded by mildly overdense gas that is still colder than the CMB and appears 
in absorption (Tg < T, STh < 0). The covering factor of material with STh < — 10 mK is 
1.7%, comparable to that of material with STf, > +10 mK: only about 1% of the pixels 
are brighter than +40 mK. While low-density gas (black color in the left-lower panel) 
remains invisible against the CMB in the NS run, the entire box glows in 21cm after being 
irradiated by X-rays. The fraction of sky emitting with STb > (+10, +20, +30, +50, +100) 
mK is now (0.57,0.31,0.19,0.1,0.035). 

The above calculations show that, even in the absence of external heating sources, 
spin exchange by H-H collisions can make filamentary structures in the IGM (heated by 
adiabatic compression or shock heating) observable in 21cm emission at redshifts z < 20. 
Some cold gas with overdensities in the range 5-100 is still detectable in absorption at a 
level of STb < — lOmK, with a signal that grows as T^^ and covers a few percent of the 
sky. X-ray radiation from miniquasars preheats the IGM to a few thousand kelvins and 
increases the electron fraction: this boosts both the H-H and the e-H collisional coupling 
between Tg and Tg, making even low-density gas visible in 21cm emission well before the 
universe is significantly reionized. Any absorption signal has disappeared, and as much 
as 30% of the sky is now shining with STi, > +20 mK. As pointed out by Nusser (2005), 
the enhanced e-H coupling makes the spin temperature very sensitive to the free-electron 
fraction: the latter is also a tracer of the H2 molecular fraction in the IGM. 



34 



Piero Madau: Early Galaxy Formation 



4.4. Concluding remarks 

Since hierarchical clustering theories provide a well-defined framework in which the his- 
tory of baryonic material can be tracked through cosmic time, probing the reionization 
epoch may then help constrain competing models for the formation of cosmic structures. 
Quite apart from uncertainties in the primordial power spectrum on small scales, how- 
ever, it is the astrophysics of baryons that makes us unable to predict when reionization 
actually occurred. Consider the following illustrative example: 

Hydrogen photo-ionization requires more than one photon above 13.6 cV per hydrogen 
atom: of order t/t 

rec ^ 10 (where tree Is the volume-averaged hydrogen recombination 
timescale) extra photons appear to be needed to keep the gas in overdense regions and 
filaments ionized against radiative recombinations (Gncdin 2000; Madau ct al. 1999). A 
'typical' stellar population produces during its lifetime about 4000 Lyman continuum 
(ionizing) photons per stellar proton. A fraction / ~ 0.25% of cosmic baryons must then 
condense into stars to supply the requisite ultraviolet flux. This estimate assumes a stan- 
dard (Salpeter) initial mass function (IMF), which determines the relative abundances 
of hot, high mass stars versus cold, low mass ones. The very first generation of stars 
('Population TIF) must have formed, however, out of unmagnctized metal-free gas: char- 
acteristics, these, which may have led to a 'top-heavy' IMF biased towards very massive 
stars (i.e. stars a few hundred times more massive than the Sun), quite different from 
the present-day Galactic case. Pop III stars emit about 10^ Lyman contimmm photons 
per stellar baryon (Bromm et al. 2001), approximately 25 times more than a standard 
stellar population. A corresponding smaller fraction of cosmic baryons would have to 
collapse then into Pop III stars to rcionizc the universe, / ^ 10^^. There are of course 
further complications. Since, at zero metallicity, mass loss through radiatively-driven 
stellar winds is expected to be negligible. Pop III stars may actually die losing only a 
small fraction of their mass. If they retain their large mass until death, stars with masses 
140 < m* < 260 Mq will encounter the electron-positron pair instability and disappear 
in a giant nuclear-powered explosion (Fryer et al. 2001), leaving no compact remnants 
and polluting the universe with the first heavy elements. In still heavier stars, however, 
oxygen and silicon burning is unable to drive an explosion, and complete collapse to a 
black hole will occur instead (Bond et al. 1984). Thin disk accretion onto a Schwarzschild 
black hole releases about 50 MeV per baryon. The conversion of a trace amount of the 
total baryonic mass into early black holes, / ~ 3 x 10^®, would then suffice to at least 
partially ionize and reheat the universe. 

The above discussion should make it clear that, despite much recent progress in our 
understanding of the formation of early cosmic structure and the high-redshift universe, 
the astrophysics of first light remains one of the missing links in galaxy formation and 
evolution studies. We are left very uncertain about the whole era from 10^ to 10" yr - 
the epoch of the first galaxies, stars, supernovae, and massive black holes. Some of the 
issues reviewed in these lectures are likely to remain a topic of lively controversy until the 
launch of the Jam,es Webb Space Telescope ( JWST). ideally suited to image the earliest 
generation of stars in the universe. If the first massive black holes form in pregalactic 
systems at very high redshifts, they will be incorporated through a series of mergers 
into larger and larger halos, sink to the center owing to dynamical friction, accrete a 
fraction of the gas in the merger remnant to become supermassive, and form binary 
systems (Volonteri et al. 2003). Their coalescence would be signaled by the emission of 
low-frequency gravitational waves detectable by the planned Laser Interferometer Space 
Antenna (LISA). An alternative way to probe the end of the dark ages and discriminate 
between different reionization histories is through 21cm tomography. Prior to the epoch 
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of full reionization, 21cm spectral features will display angular structure as well as struc- 
ture in redshift space due to inhomogeneities in the gas density field, hydrogen ionized 
fraction, and spin temperature. Radio maps will show a patchwork (both in angle and in 
frequency) of emission signals from H i zones modulated by H ii regions where no signal 
is detectable against the CMB (Ciardi & Madau 2003). The search at 21cm for the epoch 
of first light has become one of the main science drivers for the next generation of radio 
arrays. 

While many of the cosmological puzzles we have discussed can be tackled directly by 
studying distant objects, it has also become clear that many of today's "observables" 
within the Milky Way and nearby galaxies relate to events occurring at very high red- 
shifts, during and soon after the epoch of reionization. In this sense, galaxies in the Local 
Group ("near-field cosmology") can provide a crucial diagnostic link to the physical pro- 
cesses that govern structure formation and evolution in the early universe ("far-field 
cosmology" ) . It is now well established, for example, that the hierarchical mergers that 
form the halos surrounding galaxies are rather inefficient, leaving substantial amounts 
of stripped halo cores or "subhalos" orbiting within these systems (see Fig. 22). Small 
halos collapse at high redshift when the universe is very dense, so their central densities 
are correspondingly high. When these merge into larger hosts, their high densities allow 
them to resist the strong tidal forces that acts to destroy them. Gravitational interac- 
tions appear to unbind most of the mass associated with the merged progenitors, but a 
significant fraction of these small halos survives as distinct substructure. 

The Via Lactca simulation has recently shown that, in the standard CDM paradigm, 
galaxy halos should be filled with tens of thousands subhalos that appear to have no 
optically luminous counterpart: this is more than an order of magnitude larger than found 
in previous simulations. Their cumulative mass function is well-fit by A'^(> Msub) 3c ^^sub 
down to Msub = 4 X 10^ M0. Sub-substructure is apparent in all the larger satellites, 
and a few dark matter lumps are now resolved even in the solar vicinity. In Via Lactea, 
the number of dark satellites with peak circular velocities above 5 kms^^ (10 kms~^) 
exceeds 800 (120). As shown in Figure 23, such finding appears to exacerbate the so- 
called "missing satellite problem", the large mismatch between the twenty or so dwarf 
satellite galaxies observed around the Milky Way and the predicted large number of 
CDM subhalos (Moore et al. 1999; Klypin et al. 1999). Solutions involving feedback 
mechanisms that make halo substructure very inefficient in forming stars offer a possible 
way out (e.g. Bullock et al. 2000; Kravtsov et al. 2004; Moore et al. 2006). Even if 
most dark matter satellites have no optically luminous counterparts, the substructure 
population may be detectable via flux ratio anomalies in strong gravitational lenses 
(Metcalf & Madau 2001), through its effects on stellar streams (Ibata et al. 2002), or 
possibly via 7-rays from dark matter annihilation in their cores (e.g. Bergstrom et al. 
1999; Golafrancesco et al. 2006). We are coming into a new era of galaxy formation and 
evolution studies, in which fossil signatures accessible today within nearby galaxy halos 
will allow us to probe back to early epochs, and in which the basic building blocks of 
galaxies will become recognizable in the near-field. Understanding galaxy formation is 
largely about understanding the survival of substructure and baryon dissipation within 
the CDM hierarchy. 
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Figure 22. Projected dark matter density-squared map of the Via Lactea halo at the present 
epoch. The image covers an area of 800 x 600 kpc, and the projection goes through a 600 
kpc-deep cuboid containing a total of 110 million particles. The logarithmic color scale covers 
20 decades in density-square. 
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Figure 23. Left: Cumulative peak circular velocity function for all subhalos within Via Lactea's 
rvir {upper curve) and for the subpopulation within the inner O.lrvir {lower curve). Solid line 
with points: observed number of dwarf galaxy satellites around the Milky Way. Right: Same 
plotted versus virial temperature Tvir. 
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